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1.0  RESEARCH  OBJECTIVES  AND  SUMMARY 


The  work  described  in  this  report  represents  the  results  of  research  conducted  on 
turbulent  reacting  flows  in  ai r-breathi ng  propulsion  systems.  The  effort  involves 
an  assessment  of  research  on  supersonic  combustion  and  the  preparation  of  a  modular 
model  for  delivery  to  the  Air  Force  at  Wri ght-Patterson  Air  Force  Base.  The 
specific  objectives  of  this  research  are  summarized  as  follows: 

a)  SUPERSONIC  COMBUSTION 

Objective:  Under  this  task  area  an  assessment  of  supersonic  combustion 
research  was  to  be  carried  out.  The  emphasis  of  this  effort  involves  the 
identification  of  fundamental  phenomena  that  require  research  in  order  to 
improve  our  understanding  of  fuel  injection,  mixing,  flame  stabilization 
and  flame  propagation  in  supersonic  combustion  ramjet  engines. 

Research  Results:  The  results  of  the  research  conducted  on  this  task  include 
the  identification  of  pressure-turbulence  interactions  as  a  potential  expla¬ 
nation  for  trends  in  turbulent  transport  and  mixing  rates  observed  in  super¬ 
sonic  flows.  An  approach  to  the  modeling  of  these  interactions  including 
direct  coupling  of  the  turbulence  levels  to  pressure  gradient  and  pressure 
fluctuations  is  outlined  in  this  report. 

b)  MODULAR  MODEL 

Objecti ve:  Under  this  task  a  model  applicable  to  the  analysis  of  two  and 
three  dimensional  combustion  processes  is  to  be  prepared  for  the  Air  Force 
(AFWAL/PORT)  at  Wri ght-Patterson  Air  Force  Base.  A  computer  code  and 
associated  documentation  is  to  be  delivered  to  the  Air  Force. 

Research  Result:  A  basic  modular  model  computer  code  was  prepared  and 
delivered,  together  with  documentation,  to  the  Air  Force  at  Wright-Patterson 
Air  Force  Base.  This  modular  model  is  designed  for  parametric  studies  of 
the  efforts  of  operating  conditions  and  geometric  factors  on  flame  stabili¬ 
zation  and  combustion  efficiency.  It  is  specifically  relevant  to  the 
determination  of  chemical  kinetic  limitations  on  the  combustion  processes 
in  flow  regions  typically  found  in  sudden  expansion  burner  configurations  such 
as  integral  rocket  ramjet  and  ducted  rocket  combustion  chambers.  This  work 
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is  separately  documented  in  a  report  entitled  "A  Simple  Reactor 
Flow  Computer  Code  for  the  Analysis  of  Kinetic  Efforts  in  Sudden 
Expansion  Dump  Combustors",  Report  No.  SAIC 

2.0  INTRODUCTION 

The  advent  of  large  scale  space  operations  and  interest  in  wide  area  defense 
measures  have  revived  interest  in  hypersonic  propulsion  systems  using  airbreathing, 
supersonic  combustion  ramjet  engines. 

Supersonic  combustion  has  been  studied  for  a  number  of  years,  with  much  of  the  prior 
emphasis  having  been  on  the  demonstration  that  a  supersonic  fuel-air  mixture  could 
be  ignited  and  burned  under  some  conditions.  There  has  been  far  less  emphasis  on 
understanding  some  of  the  more  fundamental  physical  problems  that  accompany 
supersonic  combustion  processes.  Chemical  kinetics  and  mixing  have  been  observed 
to  be  limiting  mechanisms  in  high  speed  flows  and  have  been  responsible  for  flame 
stabilization  and  poor  combustion  problems  even  when  hydrogen  is  the  fuel.  These 
problems  become  more  critical  when  storable  fuels  such  as  liquid  hydrocarbons  and 
slurries  are  used  because  of  the  additional  times  required  for  spray  evolution, 
mixing  and  reaction  of  the  less  reactive  gaseous  and  particulate  components  of 
these  fuels.  Thus,  the  problem  is  extremely  complex  because  of  the  involvement 
of  a  number  of  interacting  mechanisms.  However,  common  to  the  supersonic  combustion 
process  and  independent  of  the  fuel  is  the  pressure  field  that  is  developed  due  to 
mixing,  as  well  as  geometry,  and  the  interaction  of  the  instantaneous  pressure  field 
with  the  turbulence  in  the  flow.  Chief  among  the  phenomena  that  are  not  well  under¬ 
stood  in  supersonic  flow  is  the  unexplained  rate  of  mixing  which  may  be  due*  at 
least  in  part,  directly  to  pressure  gradient  and  pressure  fluctuation  effects.  It 
has  been  evident  for  some  years  that  supersonic  shear  layer  mixing  rates  are 
different  from  those  encountered  in  shear  layers  at  lower  speeds,  and  indeed, 
different  from  the  mixing  rate  observed  in  a  low-speed,  two-fluid  mixii.g  layer  in 
which  the  density  ratio  is  the  same  as  that  observed  across  a  supersonic  shear 
layer.  Since  shear  layers  and  the  attendant  mixing  phenomena  are  encountered  in  a 
variety  of  forms  throughout  the  flowfield  found  within  a  supersonic-combustion 
ramjet  (scramjet)  combustion  chamber,  this  unexplained  mixing  rate  phenomenon  is 
of  critical  importance  in  understanding  the  ignition  and  combustion  process  in 
supersonic  flows. 
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Because  the  physical  cause  of  the  difference  in  mixing  rat r,  between  a  subsonic 
variable-density  (two-fluid)  mixing  layer  and  a  supersonic  mixing  layer  having 
the  same  overall  density  ratio  is  not  well  understood,  attempts  to  describe  this 
phenomenon  have  rested  primarily  on  the  development  of  empirical  factors  which 
are  intended  to  characterize  the  modification  of  the  turbulent  kinetic  energy 
distribution  in  supersonic  flows  compared  to  that  which  exists  in  an  "equivalent" 
subsonic  flow.  But  since  these  factors  are  not  founded  on  basic  grounds,  it  is 
not  surprising  that  the  application  of  them  fails  to  predict  observations  which 
have  been  made  in  supersonic  flows  of  current  interest.  For  the  same  reason, 
little  confidence  in  the  applicability  of  these  semi-empirical  factors  over  a 
broad  range  of  conditions  can  be  expressed.  Thus  what  is  required  is  the 
development  of  more  fundamental  physical  descriptions  of  the  phenomena  that  affect 
the  mixing  rate  in  supersonic  reacting  flows.  Such  a  development  will  reduce  the 
high  level  of  uncertainty  that  now  attends  estimates  of  mixing  and  combustion 
rates  in  a  scramjet  engine. 

3.0  RESEARCH  STATUS 

There  is  considerable  evidence  that  a  major  reason  for  the  lack  of  understanding 
in  supersonic  mixing  is  a  lack  of  knowledge  regarding  the  effects  of  both  mean 
and  fluctuating  pressure  gradients  on  the  mixing  rate  in  a  compressible  flow. 

These  effects  are  apparent  in  supersonic  flames  and  indeed  in  subsonic  flows, 
since  for  all  but  incompressible  flow  the  contribution  to  the  overall  turbulent 
kinetic  energy  level  of  the  pressure-related  terms  can  be  of  the  same  order  of 
magnitude  as  all  other  contributions.  Overall,  it  has  been  found  that  the 
pressure  fluctuation  and  mean  pressure  gradient  effects  can  manifest  themselves 
as  either  a  source  or  a  sink  of  turbulent  kinetic  energy.  Thus  their  effects 
on  the  overall  mixing  rate  can  be  either  positive  or  negative,  pointing  to  a 
potential  for  exercising  some  control  of  the  overall  mixing  rate  in  a  supersonic 
flow  through  control  of  pressure  fluctuations  and  mean  pressure  gradients. 

In  the  original  development  of  the  turbulent  kinetic  energy  approach  to  the 
analysis  of  turbulent  mixing  phenomena,  attention  was  restricted  to  low  speed, 
incompressible  and  non-reacting  flows.  In  this  circumstance  the  effect  of  mean 
pressure  gradients  with  respect  to  turbulence  production  is,  in  terms  of  the 
turbulent  kinetic  energy,  identically  zero.  In  an  incompressible  flow,  the  effect 
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of  a  mean  pressure  gradient  is  to  redistribute  the  fluctuating  energy  among  the 
different  components  of  the  turbulence  intensity,  without  affecting  the  overall 
energy  balance.  A  fluctuating  pressure-velocity  correlation  remains  to  be  dealt 
with,  but  since  this  has  the  effect  of  an  additional  diffusion  term,  which  is, 
in  many  cases  of  incompressible  flow,  small  compared  to  the  other  contributions 
to  the  turbulent  diffusion  of  turbulent  kinetic  energy,  this  correlation  has 
often  been  neglected.  Despite  these  rather  fundamental  limitations  with  respect 
to  the  application  of  the  turbulent  kinetic  energy  equation  to  more  complex 
compressible  and  reacting  flows,  analyses  of  these  more  complex  phenomena  have 
been  made,  often  with  unexplained  good  success.  However,  in  some  fundamental 
and  technologically  relevant  circumstances  the  analysis  of  flow  phenomena  using 
the  incompressible  turbulent  kinetic  energy  approach  has  not  been  s'.ccessful. 

The  supersonic  planar  shear  layer  provides  a  case  in  point. 

A  considerable  amount  of  data,  collected  and  analysed  for  use  in  the  1972 
NASA-Langley  conference  on  free  turbulent  shear  flows,  involved  the  planar  free 
shear  configuration.  These  data,  reviewed  by  Birch  and  Eggers  (Ref.l),  involved 
determinations  of  the  spread  rate  parameter,  (7,  in  low  speed,  single-gas,  com¬ 
pressible,  single-gas,  and  low-speed,  two-gas  flows.  This  spread  rate  parameter 
is  a  measure  of  the  width  of  the  shear  layer,  and  is  easily  obtained  from  experi¬ 
mental  velocity  profiles,  or  from  simple  visual  observations.  Although  a  is  a 
parameter  that  arises  from  the  analytical  solution  for  a  fully  self-similar  two- 
dimensional  flow,  a  specific  definition  can  be  expressed  by  the  equation 

O  =  Cl  *  (X2  -  XI)  /  (Y2  -  Yl)  (1) 


where  Y2  and  Yl  are  the  shear  layer  widths  at  X2  and  XI,  respectively.  The 
constant  Cl  depends  on  the  width  definition.  If  the  shear  layer  width,  for 
example,  is  defined  as  that  between  (U  -  U1)/(U2  -  U1 )  =  0.10  and  (U  -  U 1 ) / 

(U2  -  111 )  =  0.90,  where  U1  and  U2  are  the  velocities  on  the  different  sides  of 
the  shear  layer,  then  Cl  =  1.855  (Ref. 2).  Note  that  as  O  increases,  the  rate  of 
growth  of  the  shear  layer  decreases. 
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Birch  and  Eggers  (Ref.  1)  found  that  the  available  data  for  supersonic  shear 
layer  growth  rates  fell  into  two  distinct  groups,  as  can  be  seen  from  Fig.  1.  One 
group  of  data  indicates  a  sharp  increase  in  (J  as  the  Mach  number  increases,  while 
the  other  group,  represented  here  by  the  data  of  Hill  and  Page  (but  not,  in  fact, 
limited  to  just  these  data  alone)  shows  an  initial  rise  in  (Jas  Mach  number 
increases,  but  then  a  leveling  off  of  the  increase.  Birch  and  Eggers  suggest 
that  the  results  which  indicate  a  rapid  rise  in  the  spread  rate  parameter  with 
Mach  number  are  the  more  accurate,  since  the  results  of  Hill  and  Page  are  in  a 
developing,  rather  than  a  fully-developed  shear  layer.  However,  the  length-based 
Reynolds  number  for  these  data  is  quite  similar  to  that  of  the  other  results  shown. 
An  alternate  hypothesis  for  the  difference  in  these  results,  discussed  by  Harsha 
(Ref.  3)  notes  that  the  shear  layer  width  measured  by  Hill  and  Page  was  as  large 
as  2.4  inches  in  a  4  inch  x  4  inch  wind  tunnel,  for  the  Mach  3.5  condition. 

This  would  introduce  both  effects  of  three-dimensionality  and  of  pressure  gradients. 

The  rapid  rise  in  shear  layer  growth  rate  encountered  in  the  supersonic  shear  layer 
is  not  observed  in  the  "equivalent"  subsonic  shear  layer.  The  equivalency  referred 
to  involves  the  density  ratio  between  the  two  streams.  For  example,  if  two  air 
streams  at  the  same  total  temperature  are  considered,  one  at  M  =  2  and  the  other  at 
M  =  0.5,  a  planar  mixing  layer  between  the  streams  will  have  a  density  ratio  of  1.7 
across  it.  A  similar  density  ratio  can  be  obtained  in  a  subsonic  flow  through  the 
mixing  of  dissimilar  gases.  For  example,  the  mixing  of  streams  of  methane  and  air 
would  provide  an  equivalent  density  ratio.  But  cold-flow  mixing  of  methane  and 
air,  as  in  this  example,  would  not  result  in  the  same  shear  layer  growth  rate,  at 
the  same  velocity  and  density  ratio,  as  results  in  the  supersonic  case.  In  fact. 
Brown  and  Roshko  (Ref.  4)  found  that  a  velocity  ratio  of  0.377,  density  ratios  of 
up  to  a  factor  of  49  resulted  in  only  a  60%  change  in  (J  . 

These  observations  have  a  considerable  import  with  respect  to  the  use  of  Morkovin's 
hypothesis  in  the  application  of  turbulence  models  to  compressible  flowfield  pre¬ 
diction.  As  outlined  by  Bradshaw  (Ref.  5),  Morkovin's  hypothesis  involves  the 
argument  that  as  long  as  pressure  fluctuations  are  small  compared  to  the  mean 
pressure,  and  total  temperature  fluctuations  are  small  compared  to  the  mean  total 
temperature,  then 

P/P  — -  t'/  T~(7-  1  )  M2u'/U 


and  this  quantity  will  be  small  as  long  as  u  /  U  is  small  and  (7-  1)  M 
is  not  large  compared  to  unity.  If  the  fluctuation  in  density  is  sufficiently 
small  then  it  can  be  neglected  from  the  standpoint  of  turbulent  transport. 

Under  these  conditions  a  description  of  the  turbulent  transport  process  based 
on  incompressible  flow  is  adequate.  This  hypothesis  has  been  successfully  used 
in  the  study  of  the  development  of  compressible  turbulent  boundary  layers. 
Clearly  it  fails  for  turbulent  shear  layers,  and  because  of  the  difference 
between  supersonic  shear  layers  and  the  "equivalent"  variable-densi ty  flow,  this 
failure  is  a  result  of  flow  compressibility  alone. 


The  observation  of  compressibility  effects  in  supersonic  shear  layer  development 
led  to  the  investigation  of  several  "compressibility  correction"  techniques 
intended  to  modify  the  turbulent  kinetic  energy  distribution  in  supersonic  flow 
compared  to  that  which  exists  in  subsonic  flows.  One  such  approach  is  that 
described  by  Oh  (Ref.  6).  Oh's  approach  focused  on  the  pressure-velocity  cor¬ 
relation  term  normally  included  in  the  overall  diffusion  of  turbulent  kinetic 
energy,  and  specifically  on  a  correlation  between  the  fluctuating  pressure  and 
velocity  dilatation.  The  development  of  the  model  involves  invoking  the  exist¬ 
ence  of  eddy  Mach  waves  which,  following  Phillips  (Ref.  7)  are  supposed  to  be 
generated  whenever  the  difference  in  velocity  between  the  high-speed  side  of 
a  shear  layer  and  the  covection  velocity  of  an  eddy  exceeds  the  speed  of  sound 
in  the  high-speed  fluid.  This  model  is  sketched  in  Fig.  2,  and  requires  that 
lateral  transport  of  eddies  take  place  rapidly  enough  that  they  retain  much  of 
their  longitudinal  velocity  as  they  traverse  the  shear  layer.  Based  on  this 
model,  and  using  the  continuity  equation  to  write  the  pressure-dilatation  cor¬ 
relation  in  terms  of  mean  flow  quantities.  Oh  arrives  at  the  expression  given 
by  Eq.  2: 
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Figure  2.  Schematic  of  Eddy  Mach  Wave  Generation  Mechanism. 
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This  semi -empirical  approach  provides  the  desired  result  for  the  spread  rate  of  a 
supersonic  shear  layer,  as  shown  by  the  results  presented  in  Fig.  3;  the  restriction 
applied  to  the  use  of  the  correction  for  M<1  enables  the  subsonic  variable-density 
results  to  be  recovered  as  well.  However,  no  direction  is  provided  by  this  approach 
with  respect  to  the  treatment  of  the  supersonic,  two-fluid  shear  layer,  which  is  a 
fundamental  flow  element  in  a  scramjet  flowfield. 


Dussage,  et  al .  (Ref.  8)  consider  in  some  detail  the  properties  and  role  of  the 
production  of  turbulent  kinetic  energy  in  a  non-reacting  supersonic  flow.  In  their 
analysis,  the  turbulent  kinetic  energy  production  term  is  separated  into  three 
parts,  one  representing  production  due  to  pure  strain,  without  density  change, 
another  to  pure  dilatation  effects,  and  a  third  which  expresses  the  explicit  effects 
of  mean  pressure-gradients.  These  three  components  are  given  by  Eqs  (3)-(5)  below, 
using  mass-averaged  (Favre-averaged)  variables*: 

Pure  Strain 


-  -  I  '"K  V  1(577 


*  Favre-averaging  has  been  considered  to  provide  a  more  natural  set  of  dependent 
variables  in  a  compressible  flow.  They  can  be  defined  by  the  relation 


1  ft( 

>9  •  — 

—  Jo 


pg(t)  dt 


whereas  the  conventional  (or  Reynolds  -  average)  variable  can  be  defined  as 


a  _j_  pc 

to  Jo 


g(t)  dt 


where  g  is  any  vector  quantity  characterizing  the  flow. 


Note  that  the  fluctuating  velocity-pressure  gradient  correlation  term  in  Eq.  5  is 
neglected  in  the  rest  of  the  analysis  described  by  Dussauge,  et  al .  Using  an 
order-of-magnitude  analysis  technique,  Dussauge,  et  al .  consider  the  influence 
of  the  different  production  terms  for  a  supersonic  compression  and  a  supersonic 
expansion.  The  results  indicate  that  the  mean  pressure  gradient  production 
mechanism  enters  in  the  expanding  flow  (a  shear  layer  downstream  of  a  backstep), 
where  the  dilatation  effect  is  also  of  the  same  order  of  magnitude.  In  the  case 
of  the  compression,  the  isovolumetric  production  term  dominates  the  overall 
production  of  turbulent  kinetic  energy,  and  indeed  essentially  reduces  to  the 
"incompressible"  form.  Dussauge,  et  al .  also  consider  eddy  lifetimes  in  relation 
to  average  flow  times,  and  find  that  for  the  supersonic  expansion,  the  eddy  life¬ 
time  is  long  compared  to  the  flow  time  through  the  expansion  process.  This,  they 
argue,  indicates  that  in  this  case  the  changes  in  turbulence  pattern  that  are 
observed  must  come  from  the  production  of  turbulent  kinetic  energy  and  is  not  as 
a  result  of  alterations  in  the  turbulent  kinetic  energy  dissipation  rate,  which 
would  take  longer  to  manifest  themselves.  The  relationship  is  reversed  in  the 
compression  case  examined.  From  the  standpoint  of  its  relevance  to  a  supersonic- 
combustion  ramjet  flowfield,  this  analysis  indicates  that  different  effects  of 
mean  pressure  gradient  should  be  expected  to  be  manifest  in  different  parts  of  the 
combustor  flow. 


Substantial  effects  of  pressure  gradients  and  pressure  fluctuations  have  also  been 
predicted  to  occur  in  low  speed  flames.  For  example,  Libby  and  Bray  (Ref.  9)  show 
that  in  a  low-speed  flame  the  action  of  the  mean  pressure  gradient  or  density 
inhomogeneities  produces  countergradient  diffusion.  The  Rankine-Hugoniot  relations 
,>  show  that  the  fluid  experiences  a  small  decrease  in  pressure  upon  traversing  a 
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flame  zone.  Burned  gas  pockets,  being  of  lower  density  than  unburned  gas  pockets, 
will  respond  more  readily  to  this  decrease,  and  will  exhibit  a  tendency  to  be 
driven  farther  downstream  than  the  unburned  gas.  This  effect  on  turbulent  trans¬ 
port  can  be  dominant  (and  is  in  a  countergradient  direction)  if  the  heat  release 
is  large  compared  to  the  initial  thermal  enthalpy,  even  though  the  effect  of  the 
pressure  gradient  on  the  average  flux  of  momentum  is  entirely  negligible. 


While  Libby  and  Bray  focus  attention  on  the  effect  of  mean  pressure  gradient  and 
neglect  all  effects  of  pressure  fluctuations,  Strahle  and  Chandra  (Ref.  10) 
neglect  all  effects  of  mean  pressure  gradient  and  focus  attention  on  the  effects 
of  pressure  fluctuations,  again  for  a  low-speed  flame.  This  focus  demands  the 
development  of  modeling  of  the  pressure-velocity  correlation  term  appearing  in 
the  turbulent  kinetic  energy  balance  equation.  In  Favre-averaged  variables, 
Strahle  and  Chandra  write  the  turbulent  kinetic  energy  equation  in  the  form: 
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and  consider  in  detail  the  modeling  of  the  last  term  on  the  right-hand  side  of 
Eq.  (6).  In  equation  (6)  all  of  the  variables  are  in  a  non-dimensional  form 
which  is  itself  based  on  the  use  of  high  activation  energy  asymptotics  to  define 
the  structure  of  the  premixed  flame  being  considered.  The  term  <X>  represents  the 
dissipation  of  turbulent  kinetic  energy  and  is  neglected  in  the  further  develop¬ 
ment  of  the  modeling  in  Ref.  10.  Strahle  and  Chandra  proceed  by  dividing  the 

velocity  field  into  vortical  and  dilatational  components,  since,  loosely,  vortical 

fluctuations  can  be  identified  with  turbulence  generation  and  dilatational  fluctu¬ 
ations  with  acoustic  motions.  By  carrying  out  a  detailed  order-of-magnitude 
analysis,  they  find  that  at  high  heat  release  rates  the  pressure  fluctuations  are 
primarily  dilctational  and  hence  do  not  strongly  affect  the  turbulence  energy 
balance.  At  lower  rates  of  heat  release,  the  primary  mode  of  fluctuation  involves 

the  vortical  components,  and  it  is  under  these  circumstances  that  a  model  for  the 

pressure-veloci ty  fluctuation  term  in  Eq  (6)  is  required.  For  a  planar,  one-dimens¬ 
ional  premixed  flame,  and  under  the  somewhat  disturbing  but  nevertheless  necessary 
assumption  that  Favre-  and  Reynolds-averaged  variables  can  be  equated,  there  results 
the  expression: 
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where  a  is  a  constant  of  order  unity  and  W  represents  the  dilatational  component 
of  the  total  velocity  vector.  Because  Eq  (7)  involves  the  second  derivative  of 
the  dilatational  component  it  is  clear  that  dilatation  is  required  for  the 
correlation  to  exist.  Heat  release  is  necessary. 

By  using  Eq  (7)  in  a  suitably  modified  form  of  the  Bray-Libby  theory  (Ref.  9) 

Strahle  and  Chandra  show  that,  depending  on  the  magnitude  chosen  fora,  sub¬ 
stantial  effects  on  the  turbulent  kinetic  energy  balance  can  be  predicted  to 
arise  from  the  pressure  fluctuations  in  a  model  reacting  flow.  In  general,  Eq.  7 
represents  a  source  of  turbulent  kinetic  energy  and  acts  to  increase  the  turbulent 
flame  speed.  Taken  together,  what  both  the  analysis  of  mean  pressure  gradient 
effects  of  Ref.  9  and  of  pressure  gradient  fluctuations  in  Ref.  10  show  is  that 
neglecting  either  of  these  terms  in  a  general  reacting  flow  problem  can  be  an 
erroneous  assumption.  However,  both  of  these  analyses  are  limited  to  a  simple 
one-dimensional  premixed,  fast  chemistry,  flame  approximation,  and  a  great  deal 
of  work  remains  to  be  done  to  generalize  these  results  to  conditions  of  interest 
in  the  analysis  of  supersonic  reacting  flow. 

An  important  factor  in  a  theoretical  approach  is  the  role  played  by  the  turbulent 
length  scale.  Representing  a  mean  eddy  size,  this  length  scale  enters  the 
definition  of  the  turbulent  kinetic  energy  dissipation  rate,  and  has  been  utilized 
as  well  in  the  definition  of  modeling  for  the  fluctuating  pressure-strain  cor¬ 
relation  terms  (e.g.,  Rotta,  Ref.  13).  However,  in  all  turbulent  shear  flows,  a 
spectrum  of  eddy  sizes  and  therefore  length  scales  exists,  and  thus  the  use  of  a 
single  length  scale  is  equivalent  to  an  assumption  that  all  length  scales  are 
proportional,  so  that  any  one  can  be  used  to  characterize  the  flow  -  an  "equilibrium 
assumption.  Our  hypothesis  is  that  in  the  presence  of  strong  pressure  gradients 
and  compressibility  effects  this  "equilibrium"  no  longer  exists,  so  that  different 
length  scales  need  to  be  defined  for  turbulence  dissipation  and  for  the  pressure- 
strain  interaction,  for  example. 

One  piece  of  evidence  for  this  hypothesis  is  the  work  reported  by  Hanjalic',  et  al . 
(Ref.  14)  involving  the  passage  of  ?  grid-generated  decaying  turbulent  flow  through 
a  contraction  section.  During  this  process,  the  kinetic  energy  is  observed  to 
increase  and  then  decay  again,  but  at  a  slower  rate  than  before.  Two-equation 


model  predictions  either  underpredict  the  secondary  energy  peak  and  the  subsequent 
rate  of  decay,  or  overpredict  both,  depending  on  the  constants  selected.  A  model 
incorporating  two  scales  does  successfully  predict  this  flowfield,  and  the  results 
show  that  as  the  flow  proceeds  through  the  contraction  section  (and  thus  through 
a  mean  pressure  gradient),  the  initially  equilibrium  relationship  between  the  two 
scales  of  motion  is  altered. 

An  investigation  was  carried  out  of  the  multiple-dissipation-length  scale  approach 
(Fabris,  Harsha,  and  Edelman)  Ref.  15  and  it  found  that  even  in  relatively  simple 
flows  it  provides  a  greater  generality  than  the  basic  two-equation  model:  round 
jet  predictions  can  be  made  with  the  same  coefficients  as  for  other  simple  flows, 
for  example.  However,  this  work  did  not  involve  consideration  of  the  pressure- 
strain  term,  or  of  strong  pressure  gradients,  or  of  compressible  flow,  all  of 
which  are  of  potential  importance  in  high  speed  flows. 

Viewed  from  the  perspective  of  a  multiple  length  scale  model,  the  existence  of 
coherent  structures  represents  the  behavior  of  the  flow  at  the  largest  scales  of 
interest.  Energy  is  transferred  from  these  scales  to  the  smaller  scales  of 
motion,  and  it  is  this  transfer  process  that  provides  the  connection  between  the 
observed  coherent  structure  and  classical  turbulent  flow  models.  Previous  work 
(Hanjalic',  et  al .  Ref.  14  ,  Fabris,  et  al.  Ref.  15)  has  limited  the  number  of 
scales  to  be  considered,  the  transfer  mechanism  which  we  believe  is  Reynolds 
number  dependent  and  the  equation  formulation.  Thus,  in  the  previous  work  the 
Reynolds  averaged  variables  have  been  used,  with  two  scales  and  a  direct  transfer 
of  energy  independent  of  local  flow  Reynolds  number.  For  the  flows  to  be  considered 
in  this  work,  the  approach  is  to  systematically  increase  the  number  of  scales 
considered,  with  the  largest  scale  described  as  a  periodic  (or  quasi-periodic) 
fluctuation,  and  a  local  Reynolds  number-related  energy  transfer  rate.  At  the 
other  end  of  the  spectrum,  a  small-scale  energy  dissipation  rate  can  be  defined 
to  represent  mixing  on  a  scale  consistent  with  that  which  applies  to  chemical 
reactions . 

In  contrast,  both  the  algebraic  Reynolds  stress  model  (ASM)  and  the  k-e  approaches 
are  single-point  models  which  adopt  a  single  time  scale  proportional  to  the  turb¬ 
ulence  energy  turnover  time,  k/e  .  However,  it  is  overly  simplistic,  at  least 
conceptually,  to  assume  that  a  single  time  scale  can  successfully  characterize  the 
rates  of  progress  of  different  turbulent  interactions. 


The  key  to  the  new  multiple-scale  approach  is  the  recognition  that  while  the 
dissipation  equation  and  the  kinetic  energy  equation  both  contain  production 
and  dissipation  terms,  these  processes  occur  in  different  spectral  regions  of 
the  flow.  That  is,  turbulence  energy  production  occurs  in  the  larger  eddies 
in  the  flow,  while  dissipation  phenomena  involve  primarily  the  smaller  scales. 

Thus,  there  must  be  a  transfer  of  energy  from  the  larger  scales  to  the  smaller, 
and  this  transfer  can,  in  certain  situations,  introduce  a  lag  phenomenon,  so 
that  turbulence  energy  production  and  turbulence  dissipation  do  not  necessarily 
both  increase  or  decrease  in  the  same  region  of  the  flow  as  is  implied  by  a 
single-scale  model . 

To  introduce  a  model  in  which  the  evolution  of  the  different  scales  appropriate 
to  the  large-eddy  production  region  and  the  small-eddy  dissipation  region  can  be 
accounted  for,  we  introduce  a  partitioning  of  the  turbulence  energy  and  its 
dissipation  rate  into  three  regions.  For  wave  numbers  less  than  k^,  a  production 
region  is  defined,  characterized  by  a  turbulent  kinetic  energy  kp  and  a  dissipation 
rate  €p.  This  dissipation  rate  controls  the  transfer  of  energy  through  the  trans¬ 
fer  region  k^c/cck For  wave  numbers  higher  than  turbulence  energy  is 
dissipated  as  heat.  A  separate  kinetic  energy  and  dissipation  rate  equation  is 
written  for  the  transfer  region,  characterized  by  k^  and  e T,  and  the  production 
term  in  the  kinetic  energy  equation  for  the  transfer  region  is  equal  to  the 
dissipation  rate  ep  in  the  production  region. 

\ 

The  partitioning  of  the  energy  spectrum  into  three  regions,  as  in  this  example 
requires  two  sets  of  transport  equations,  given  the  assumption  (basic  to  most 
turbulence  modeling)  that  the  mechanisms  involved  in  the  final  dissipation  of 
turbulent  kinetic  energy  into  thermal  energy  are  capable  of  accepting  all  of  the 
energy  transferred  to  them.  This  assumption  is  the  reason  that  the  physical  fluid 
viscosity  does  not  appear  in  the  turbulence  dissipation  rate  equations. 

The  model  equations  for  the  production  and  transfer  region  turbulent  kinetic 
energy  and  dissipation  rate  are  similar  in  form  to  the  standard  k  and  €  transport 
equations.  They  can  be  written  as  follows: 
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in  which  the  subscript  p  refers  to  the  production  region  and  T  to  the  transfer 
region.  In  this  formulation,  the  turbulent  viscosity  is  given  by 


In  this  example,  the  equations  are  written  in  incompressible  Reynol ds-averaged 
form.  For  a  compressible  flow,  these  equations  can  be  derived  either  in  Reynolds- 
averaged  or  Favre-averaged  form.  If  they  are  derived  in  Reynolds-averaged  variables, 
correlations  which  involve  density  fluctuations  are  introduced,  and  these  turbulence 
correlation  expressions  take  a  different  form  if  Favre-averaged  variables  are  used. 
The  choice  of  which  form  of  equations  to  use  in  the  theoretical  analysis  is  depend¬ 
ent  on  the  instrumentation  and  measurements  involved  in  the  experimental  program: 
laser  diagnostics,  for  example,  measure  the  time-average  fluctuations  which  can  be 
used  to  construct  the  Reynol ds-average  correlations,  providing  a  means  for  obtaining 
density  fluctuations  and  their  correlations  (as  well  as  pressure  fluctuations)  is 
available.  In  addition,  the  equations  have  been  written  in  time-dependent  form  to 
account  for  large-scale,  quasi-periodic  motions  of  coherent  structures. 

Note  also  that  in  writing  equations  8-11  the  pressure  fluctuation  terms  have  been 
neglected,  but  their  inclusion  is  straight  forward.  Also,  equations  8-11  indicate 
a  two-scale  breakdown  in  which  the  production  of  transfer  region  kinetic  energy  is 
directly  the  dissipation  of  production  region  kinetic  energy.  In  actuality,  there 
are  many  more  steps  to  the  energy  cascade  process.  For  example,  energy  spectrum 
observations  indicate  that  the  wave  number  range  over  which  significant  energy  is 
contained  in  the  turbulence  is  considerably  wider  in  the  interior  region  of  the 
flow  than  on  the  edges,  suggesting,  at  the  outer  edge  of  the  flow  which  is  dominated 
by  the  large  eddies,  a  more  direct  transfer  of  kinetic  energy  from  the  large  scales 
to  the  dissipating  scales  than  is  the  case  in  the  flowfield  interior.  Since  it  is 
also  observed  that  the  large  structure  is  more  evident  at  low  Reynolds  numbers  than 
at  higher  Reynolds  numbers  in  jet  and  shear  layer  flows,  it.  is  also  clear  that  the 
energy  transfer  mechanism  is  Reynolds-number  dependent. 

In  order  to  gain  insight  on  the  magnitude  and  coupling  of  the  pressure  field  in 
supersonic  ducted  flow  a  baseline  analysis  of  a  simple  supersonic  combustor  flow 
field  has  been  carried  out. 

Common  to  the  variety  of  engines  used  for  hypersonic  flight,  such  as  rocket  engines, 
airbreathing  engines  and  composite  engines  (combinations  of  both)  are  the  problems 
associated  with  propellant  combustion.  With  the  exception  of  the  pure  solid 
rocket,  the  combustion  process  is  coupled  to  a  fuel /oxi di zer  mixing  process.  If 
a  spectrum  of  pressures,  temperatures  and  flow  velocities  are  considered,  the 
combustion  process  may  be  (1)  diffusion  controlled,  i.e.,  fuel  dependent  only 


upon  the  mixing  rate;  (2)  reaction  controlled  wherein  the  kinetics  of  the 
oxidation  process  determine  the  heat  release  rate;  and  (3)  an  intermediate 
regime  wherein  the  heat  release  rate  depends  upon  the  coupled  mixing  and  kinetic 
rates . 

The  diffusion  controlled  process  is  desirable  in  comparison  to  the  reaction 
controlled  situation  because  the  latter  involves  mixing  process  followed  by  the 
combustion  process  and  could  lead  to  excessively  long  combustion  chambers.  In 
addition,  for  hypersonic  flight  involving  pure  airbreathing  and  composite  engines, 
shock  losses,  dissociation  losses  and  heat  transfer  limitations  require  that  the 
flow  remain  supersonic  during  the  combustion  process.  Analysis  of  the  supersonic 
combustion  process  requires  consideration  of  the  coupling  between  the  mixing  and 
kinetic  processes  because  flow,  mixing  and  reaction  times  are  generally  of  equal 
order  of  magnitude.  Of  particular  importance  in  high  speed  combustion  chamber 
flows  is  the  optimization  of  the  combustion  chamber  length  required  to  achieve 
efficient  heat  release  while  maintaining  reasonable  chamber  lengths.  To  achieve 
an  understanding  of  the  relevant  parameters  and  their  inter-relationships,  it  is 
necessary  to  analyze  the  details  of  the  reacting  flowfield.  A  meaningful  analysis 
for  such  flow  configurations  requires  including  turbulent  mixing,  finite  rate 
combustion  and  pressure  gradients  in  a  model  which  appropri ately  accounts  for  the 
coupling  of  these  phenomena.  It  should  be  noted  that  the  pressure  field  must  be 
considered  in  detail  because  of  the  sensitivity  to  pressure  of  both  the  chemical 
kinetics  and  velocity  field.  The  structure  of  the  pressure  field  includes  the 
locally  induced  pressure  due  to  mixing  and  combustion  as  well  as  the  pressure 
impressed  on  the  flow  due  to  the  interactions  with  the  confining  walls. 

A  configuration  which  characterizes  the  flows  considered  in  this  investigation 
is  shown  in  Fig.  4.  Although  the  complete  flowfield  can  involve  a  complex  internal 
shock  structure  the  present  form  of  the  analysis  does  not  resolve  the  pressure 
field  in  terms  of  discrete  shock  structures. 

The  classical  analysis  of  viscous  flow  problems  involving  unequal  characteristic 
length  scales  in  the  coordinate  directions  leads  to  the  boundary-layer  formulation. 
The  describing  equations  are  parabolic  and  the  pressure  is  constant  across  the 
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smaller,  or  lateral,  dimension  of  the  flow.  The  pressure  variation  in  the  primary 
flow  direction  is  either  specified  as  a  given  distribution  or  is  determined  by 
approximate  methods  based  on  the  induced  pressure  due  to  the  displacement  of  the 
inviscid  outer  flow.  In  either  case  the  pressure  across  the  flow  is  assumed 
constant.  In  flows  involving  combustion  an  accurate  prediction  of  the  induced 
pressure  field  is  crucial  to  predicting  the  force  on  solid  surfaces  adjacent  to  the 
combustion  region.  Moreover,  the  ignition  delay  and  reaction  time  of  the  combustion 
process  itself  may  be  critically  dependent  on  the  local  pressure,  particularly  for 
pressure  levels  below  atmospheric.  Furthermore,  non-uniform  transverse  pressure 


distributions  are  inherent  in  jets  operating  either  in  an  under-expanded  or  over- 
expanded  configuration.  Finally,  the  coupling  of  the  pressure  field  to  the 
turbulence  levels  included  in  the  supersonic  flow  requires  that  the  pressure  field 
be  accurately  predicted. 

The  analysis  of  this  problem  requires  a  study  of  the  full  Navier-Stokes  equations. 
The  approach  involves  the  choice  of  coordinates  and  a  set  of  physically  reasonable 
assumptions  which  will  provide  a  tractable  formulation  and  yet  preserve  the  basic 
character  of  the  flow. 

The  general  describing  equations  for  the  steady  flow  of  a  reacting  mixture  of 
perfect  gases  are  given  by: 

Global  Continuity: 

V  •  p V  =  0  (14) 

Species  Continuity: 

V  •  P« ,-V  =  P  -  V  •  j1  (15) 

i  =  ith  specie 

Where  j^  is  the  diffusional  mass  flux  vector: 
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Momentum: 


V-pVV  =  -  P+V-r 


Where  t  is  the  viscous  stress  tensor. 


:nerqy: 


y.pVH  =  V  •  (  T  •  V)  -  V  •  3  -  V  •  ?h  .j.  (17) 

V  .  V 

Where  H  =  ?<Xh.  (T)  +  -  is  the  stagnation  enthalpy,  and  q  is  the  heat 

ill  2 

conduction  vector. 

State: 
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Transport  Processes 

The  transport  processes  considered  in  this  investigation  include  viscous  stresses, 
species  diffusion  and  heat  conduction.  The  models  used  to  represent  these  processes 
are  as  follows: 

Vicous  Stress 

T  =  (  -  2/3  5  :  e  6  (19) 


Where 


e  =  1/ 2 ( vv  +  v  *  v) 


and  6  is  the  unit  tensor. 


Species  Diffusion 

The  diffusion  of  each  species  is  assumed  to  depend  only  on  the  gradient  of  the 
particular  species  mass  fraction.  This  assumption  requires  equal  binary  diffusion 
coefficients  D. .,  for  each  species  so  that  Fick's  law  is  applicable. 

'  J 


Under  this  assumption  the  diffusional  mass  flux  is  given  by: 
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P°no 


Vc^  = 
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(21) 


Heat  Conduction 

Consistent  with  the  above  assumptions  the  heat  conduction  may  be  represented 

^-PVT  (22) 


A  schematic  of  the  coordinate  system  adopted  in  this  analysis  is  shown  in 
Fig.  5.  A  coordinate  system  is  sought  wherein  a  minimum  of  physical  assumptions 
are  required  to  simplify  the  Navier-Stokes  equations  and  within  which  a  numerical 
solution  is  most  easily  formulated  and  implemented.  The  implicit  coordinate 
system  provides  these  desired  conditions.  Some  insight  into  this  choice  is 
obtained  by  considering  the  important  case  of  an  underexpanded  (or  overexpanded) 
jet.  The  initial  flow  deflection  is  governed  almost  entirely  by  the  adjustment 
of  the  pressure.  The  jet  flow  is  separated  from  the  external  flow  by  the  "dividing" 
streamline  and  the  gradients  in  velocity,  concentration  and  temperature  will  be 
greatest  in  the  direction  normal  to  the  streamline.  Of  course,  as  the  flow  becomes 
more  parallel  to  the  jet  axis  this  condition  is  still  satisfied.  It  appears  then 
that  in  this  coordinate  system  some  of  the  usual  boundary-layer  assumptions  may 
be  applied.  Furthermore,  it  is  shown  later  that  a  natural  consequence  of  this  is 
a  characteristic  method  of  solution  which  can  be  implemented  with  a  certain  degree 
of  ease. 

Edelman  et  al..  Ref.  16  (1968)  have  made  a  systematic  order  of  magnitude  analysis 
in  the  s-n  coordinate  system  (direction  along  and  normal  to  a  streamline)  for  Eqs. 
(14)  to  (22)  for  both  two-dimensional  and  axisymmetric  flows.  The  problem  involves 
a  class  of  supersonic  jet  flows  having  a  principal  flow  direction  which  are 
conditions  also  found  in  practical  configurations  including  rocket  nozzles  and 
supersonic  combustion  ramjet  engine  combustion  chambers.  In  a  supersonic  flow, 
small  deflections  can  induce  a  non-uniform  pressure  field  which  can  have  a 


by  the  Fourier  law  given  by: 


a  = 


k  VT 


where 


T,a 

i  i 


significant  effect  upon  combustion  and  pressure  forces  on  surfaces  adjacent  to  the 
flow.  Furthermore,  if  the  flow  curvature  is  small  in  regions  where  viscous  effects 
dominate  across  the  entire  field,  viz.,  away  from  the  near  region  of  an  under¬ 
expanded  jet,  then  boundary  layer  type  approximations  are  applicable  to  the  diffusiv 
terms  in  the  describing  equations.  They  have  shown  that  Eqs.  (14)  to  (22)  can  be 
reduced  to  the  following  simple  form. 


Global  Continuity: 


(pq)s  ♦  \Pq'  sin  0  +  (pq)n=  0 


s-Momentum: 
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n-Momentum: 


PP  es  +  Pn  “  0 


Species  Continuity: 


Pi  <ai>s  ‘  p“i +  7J 


J„Le  (a,), 


Energy: 
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In  equations  (23)  -  (29),  the  subscript  notation  has  been  used  to  indicate 
partial  derivatives,  i.e., 

(pq)s  -~k  <pq> 


rrj  JL  H  1  =  _L.rri£_  _SH] 
Lr  Pr  nJ  n  3n  [rPr  ^n  J 


The  solution  technique  applied  to  equations  23  through  27  involves  the  coupling 
of  both  hyperbolic  and  parabolic  equations.  It  is  the  purpose  of  the  following 
discussions  to  describe  the  details  of  the  method  of  solution  which  can  be 
described  in  terms  of  an  "inner"  calculation  (characteristics)  and  an  "outer" 
calculation  (mixing). 

The  characteristics  calculation  is  initiated  in  a  standard  manner  from  a  "line" 
of  given  data  as  shown  in  Fig.  6. 

The  characteristic  directions  are  dependent  only  upon  the  local  state  and  the  left 
(I)  and  right  (II)  running  characteristic  directions  from  point  K,  say,  are  given 
by: 


Tan  eK  -  + 


-Vp-TiTTT 


Where 


F  =  T 


2  a 


dh.  (T) 


and  all  properties  are  evaluated  at  the  point  K. 
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The  compatibility  equations  are  obtained  in  the  usual  way  assuming,  however,  that 


the  viscous  terms  are  part  of  the  forcing  function.  These  equations  provide  a 

_  .  .  i  J..  .  ,  •  /a  T  1  /r  \ J  1.  I _  .  . _  n  .  .  m 


relation  between  the  flow  deflection  (0=  Tan  T  )  and  the  pressure  p,  and  are 
given  by: 


dT  +  .  dp  + 


Tf-  1 


Where 


l  +  t‘ 


pq  Tanf 


(1 +r  )sine 


F  ,  /  .  q„  C0S» 


2  h  D  + 

1  1  1  P 


-JL  a. 
pq  inn 


+  _c_ 

pq 


P_  cosfl  a  +  _i_ 
>q  y  in  q 


The  chemical  production  term,  ,  represents  the  production  rate  of  species  i 
through  the  mechanism  of  a  system  of  chemical  reactions  given  by: 


s  rjflj 
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m  =  1 ,  M 


where  N  species  are  involved  in  M  reactions  and  the  p.'s  are  the  stoichiometric 
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coefficients  associated  with  Aj  species.  The  production  rate,  ,  is  then  given  by: 
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is  the  molar  density  of  the  i  specie. 

The  sequence  of  calculations  is  initiated  with  a  characteristic  type  calculation 
Referring  to  Fig.  6,  the  up-running  characteristic  from  point  K-l  and  the  down¬ 
running  characteristics  from  point  K  are  constructed  and  their  point  of  inter¬ 
section,  CK>  is  obtained.  This  procedure  is  carried  out  for  all  intersections 
across  the  field.  The  minimum  forward  distance  required  for  an  intersection 
determines  the  "characteristics"  step  size  which  guarantees  that  all  new  data 
computed  across  the  field  is  within  the  domain  of  influence  of  the  initial  data. 

The  next  step  in  the  sequence  of  calculations  involves  the  solution  of  the 
compatibility  equations  along  the  characteristics.  The  differenced  form  of  Eqs. 
(32)  are  given  by: 
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Solving  Eqs.  (39)  and  (40)  for  the  pressure  gives: 

Vl  -TK  +  9K-1  +  9K^II  +  bK-lPK-l  +  bK  PK 


(43) 


bK-l  +  bK 


Tr  is  determined  by  substitution  into  either  Eqs.  (39)  or  (40) 


As  discussed  above,  a  single  forward  step  size  is  used  and  this  is  done  in 
anticipation  of  the  subsequent  mixing  calculation  discussed  below.  However, 
this  requires  interpolations  for  the  data  required  at  points  A,  D  and  finally 
E  (Fig.  6).  Thus,  the  pressure,  p,  and  flow  deflection,  T,  are  determined  at 
A  and  D  by  interpolation  along  the  characteristics  and  then  a  second  interpolation 
is  performed  across  the  field  to  determine  p  and  T  on  the  streamline  at  point  E. 

In  this  way,  the  pressure  gradient  along  the  streamline  is  obtained: 
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(44) 


This  result  is  utilized  in  the  mixing  calculation  discussed  in  the  next  section. 


The  explicit  bilateral  coupling  between  the  pressure  field  and  the  viscous  mixing 
process  is  through  the  forcing  function  g  in  Eqs.  (39)  and  (40)  and  the  momentum 
equation  along  the  streamline,  respectively.  The  former  contains  the  diffusion 
and  chemical  reaction  terms  and  the  latter  contains  the  pressure  gradient. 
Furthermore,  the  description  given  above  for  the  p-T  determination  represents  a 
single  calculation  in  an  iterative  calculation.  Thus,  the  forcing  function  is 
modified  by  averaging  it  over  a  step  which  includes  the  "new"  diffusion  contributions 
obtained  at  the  end  of  the  mixing  calculation.  The  p-T  calculation  is  repeated 
until  the  difference  in  two  successively  determined  pressures  is  less  than  a 
prescribed  tolerance. 


The  calculation  procedure  as  described  so  far  lends  itself  uniquely  to  the  appli¬ 
cation  of  already  existing  numerical  solutions  of  the  parabolic  mixing  problem. 
The  solution  is  obtained  by  an  explicit  finite  difference  representation  of  the 
describing  equations  along  the  streamlines.  These  equations  are  given  by: 
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It  is  important  to  note  that  the  describing  equations  for  mean  quantities  have 
been  simplified  from  their  most  general  form  by  the  assumption  that  the  Prandtl 
and  Schmidt  numbers  are  unity.  Furthermore,  in  the  numerical  representation  of 
the  above  equations  it  is  assumed  that  n  — y.  Both  of  these  approximations  were 
made  to  facilitate  the  initial  stages  of  the  program  development  but  it  is  believed 
that  the  limitations  are  not  critical. 


To  illustrate  the  numerical  solution  technique,  consider  the  difference  form  of 
Eq.  (45)  with  reference  to  Fig.  7. 
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where 


is  given  by  the  characteristics  calculation  (see  Eq.  (44). 


Eq . (51 )  represents  an  explicit  finite  difference  form  of  the  s-momentum  equation. 

The  energy  and  species  conservation  equations  assume  a  similar  form.  The  forward 
step  size  is  limited  by  a  stability  requirement  which  for  linear  parabolic  equations 
is  given  by: 


As  <  1/2  pq 


To  insure  stability  of  the  non-linear  system  a  fraction  of  the  characteristic  step, 
O,  is  used: 


where  Ns  is  an  input  integer.  In  cases  involving  fast  reacting  systems  Ns=6  has 


been  found  to  be  sufficient. 


The  species  conservation  equation,  Eq.  (47),  is  treated  in  a  special  way  for  configu¬ 
rations  involving  finite  rate  kinetics.  The  equations  for  turbulent  kinetic  energy 
e  and  dissipation  rate  e  are  also  treated  differently  to  avoid  negative  values 
which  are  unphysical.  In  the  beginning  of  the  calculations  when  both  diffusion 
and  production  are  negligible,  the  dissipation  term  (last  terms  in  Eqs.  (48)  and 
(50)  when  expressed  explicitly  in  the  finite  difference  representation  can  give 


negative  values  for  e  and  e  under  certain  circumstances.  This  is  avoided  by  partly 
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implicitizing  the  dissipation  term.  Thus,  e  and  e  in  the  dissipation  term  are 


approximated  as  (em  „,,)  (em  )  u  and  (  e  )  (  e  )  respectively.  Thus  Eq. 

m,n+l'  m,n  m,n+l  m,n  K  J  M 


(48)  can  be  approximated  as: 


em,n+l  =  em,n  +  ^m,n  /  1+ 


Vs  a0  e*  „ 
_ _ 2  m,n 


q  £  . 
vn,n  *  k 


(54) 


where  X  represents  the  diffusion  and  production  term. 
m,n 

The  above  discussions  pertain  specifically  to  the  solution  technique  for  the 
describing  partial  differential  equations.  To  completely  specify  a  given  problem 
the  initial  and  boundary  conditions  must  be  given.  In  the  present  report  a  shock¬ 
less  ducted  jet  flow  configuration  is  described.  The  wall  is  assumed  to  be 
impermeable,  and  adiabatic.  The  boundary  conditions  are  given  by: 
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The  arbitrary  wall  contour,  y  =  y  (x)  may  be  specified  in  terms  of  the  coefficients 

t  h  ^ 

of  5in  order  polynominals  covering  four  adjoining  regions: 


y  =  X  a  x 
w  0  n 


n 


(56) 


The  pressure  is  determined  from  the  appropriate  uprunning  characteristic  since 
the  flow  deflection  at  the  "new"  wall  point  is  known  from  Eq.  (56). 


The  calculation  proceeds  as  follows  (see  Fig.  8) 


A  characteristic  step  is  taken  to  determine  the  minimum  forward  distance,  O  ,  for 

intersections.  A  star,  *,  point  is  found  by  interpolation  between  the  w  and  KT 

n  n 

points  whose  uprunning  characteristic  intersects  the  wall  at  the  end  of  the  step, 
xn+i=  xn  +  a.  The  pressure  at  the  wall  is  then  computed  through  the  compatibility 
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equation: 
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The  remaining  properties  along  the  n+1  line  are  determined  and  the  diffusion 
calculation  is  performed.  The  wall  boundary  conditions  for  the  diffusion  step, 
Eqs.  (55)  are  imposed  by  inserting  a  data  point  above  the  wall  denoted  by  KT+1. 
The  numerical  difference  form  of  Eqs.  (55)  are  then  given  by: 


a.  =  a. 

1 KT+1  7KT 


hkt+i 

hkt 

(58) 

—  n  _  _ 

Cf 

(pq)">ean  /  \ 

qKT+l 

qKT 

2 

H  ( y KT+1  yKT ) 

where  the  physical  location 

of  the  KT+1 
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is  given  by: 

yKT+l  =  yw  ‘  yKT 


Finally,  the  iterative  process  between  the  characteristic  step  and  the  diffusion 
step  is  performed  as  described  before. 

This  analysis  has  been  applied  to  a  jet  flow  configuration  bounded  by  a  parallel 

wall.  This  configuration  is  basic  to  many  fuel  injection  problems.  The  results 

for  an  underexpanded  all-supersonic  ducted  jet  with  P.fP  =2  and  with  hydrogen 

J  — 

(primary  flow)  exhausting  in  air  (secondary  flow)  are  presented  in  this  section; 
initial  conditions  for  this  calculation  are  given  in  Table  1.  Two  sets  of 
calculations  have  been  carried  out;  one  set  with  a  constant  eddy  viscosity  and 
the  second  set  with  a  viscosity  computed  by  the  two-equation  TKE  model.  Figs.  9a 
through  9f  show  the  results  of  the  first  case  where  /t  =  0.00449  lb-sec/ft  was 
chosen,  ihis  is  a  rather  high  value  for  the  viscosity  and  the  results  reflect 
this  fact.  Figure  9a  shows  the  wall  and  axis  pressure  distribution  up  to 
x/r,=  44.  The  initial  rise  in  pressure  on  the  axis  is  related  to  the  value  of 


TABLE  1. 


INITIAL  CONDITIONS  FOR  SAMPLE  CASE 


AXISYMMETRIC  HYDROGEN  INJECTION  INTO  AIR  STREAM 


Stream 

Mach  No. 

Veloci ty 

Pressure 

Temperature 

(Ft/sec) 

(LBF/Ft2) 

(°k) 

Air 

3.375 

8333 

5290 

1500 

Hydrogen 


2.032 


16,634 


10,580 


1100 


9  b.  Axial  Decay  of  Hydrogen  Mass  Fraction. 
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viscosity  and  the  initial  pressure  ratio.  Later  with  the  adoption  of  the  two- 
equation  TKE  model  this  sharp  rise  disappears.  The  axial  decay  of  hydrogen  mass 
fraction  is  shown  in  Fig.  9b.  Because  of  the  high  viscosity  the  decay  is  almost 
immediate  with  no  potential  core  formation.  Again,  later  we  shall  show  that  this 
shortcoming  is  eliminated  by  the  use  of  the  TKE  model  for  turbulence  prescription. 
Figures  9c  to  9f  show  the  profiles  of  pressure,  velocity,  temperature  and  hydrogen 


mass  fraction  at  x/r. 

J 


40.6,  which  are  typical  of  a  well  mixed  flowfield. 


Figs.  10a  to  lOf  represent  the  flowfield  for  similar  initial  conditions  and  geometry 
but  with  the  constant  viscosity  replaced  by  the  two-equation  TKE  model.  Fig.  10a 
shows  the  distribution  of  wall  and  axis  pressure.  The  initial  sharp  rise  in  the 
centerline  pressure  present  in  the  previous  case  goes  away  with  the  adoption  of 
this  turbulence  model.  Fig.  10b  shows  the  axial  decay  of  hydrogen  mass  fraction. 
Based  on  this  the  potential  core  length  can  be  estimated  to  be  approximately  20 
jet  diameters.  Figs.  10c  to  lOf  show  the  profiles  of  pressure,  velocity,  temper¬ 
ature  and  hydrogen  mass  fraction  at  x/r.=  38.4,  which  in  general  show  that  the 

J 

width  of  the  pressure  field  is  substantially  greater  than  the  spread  of  the  viscous 
region. 


i 
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The  results  of  these  calculations  show  the  possibility  of  the  development  of 
relatively  large  pressure  gradients  in  supersonic  combustor  flows.  The  sensitivity 
of  the  induced  pressure  field  to  turbulence  model  assumptions  is  also  demonstrated 
by  these  results.  However,  conventional  turbulence  modeling  was  used  in  this 
analysis  and  the  direct  effect  of  the  pressure  gradient  on  turbulence  level  was  not 
included  in  these  calculations.  Furthermore,  the  explicit  effects  of  pressure 
fluctuations  on  the  turbulence  was  also  neglected  in  these  baseline  calculations. 
These  pressure  coupling  effects  are  the  key  elements  that  have  been  previously 
cited  in  this  report  as  requiring  inclusion  in  the  analysis  and  the  results  of  the 

Preliminary  calculations  indicate  that  the  induced  pressure  gradients  can  play  a 
significant  role  in  the  turbulence-pressure  interaction  phenomenon. 


Because  reacting  flow  is  the  ultimate  application  of  the  research  to  be  carried 
>'  out  in  this  program,  the  chemical  kinetics  of  the  reaction  process  must  also  be 
considered  in  the  analysis.  Short  residence  times  and  relatively  low  pressures 

m  •« 

y  such  as  are  encountered  in  a  scramjet,  increase  the  importance  of  chemical  kinetics 
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Two  Equation  TKE  Model 


Axial  Decay 
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Figure  10  c.  Pressure  Profile  at  x/R  =  38.4.  Figure  10  d.  Velocity  Profile  at  x/R 


Two  Equation  TKE  Model 


Figure  10c.  Temperature  Profile  at  x/R.=38.4.  Figure  lOd.  Hydrogen  Mass  Fraction  Profile  at 

J  x/R  .=38.4 . 


effects  in  a  scramjet  relative  to  what  is  observed  in  subsonic  combustors.  At 
the  same  time,  turbulence-chemistry  effects  also  become  important.  The  analysis 
of  the  interaction  between  turbulence  and  chemical  reaction  rates  is  a  subject 
of  considerable  current  interest,  but  no  tractable  approach  for  including  a 
general  and  detailed  finite-rate  kinetics  model  in  the  analysis  of  turbulence- 
combustor  interactions  has  yet  been  developed.  Although  the  relatively  small 
pressure  fluctuations  and  pressure  gradients  strongly  affect  the  turbulent 
kinetic  energy  balance  their  influence  on  chemical  kinetic  rates  needs  to  be 
determined.  Thus,  as  a  first  step  in  the  consideration  of  all  of  the  effects 
involved  in  turbulent  mixing  and  reaction  in  a  supersonic  flow,  a  fast-chemistry, 
probability  distribution  function  closure  for  the  reaction  problem  should  be 
used.  However,  stronger  pressure  gradients  can  be  used  to  control  the  combustion 
process  in  a  supersonic  flow,  and  these  will  also  affect  the  mixing  process.  Thus 
an  assessment  of  the  effects  of  pressure  gradients  and  pressure  fluctuations  in  an 
environment  where  chemical  kinetics  effects  are  important  must  also  be  considered. 


5.0  INTERACTIONS 


The  sudden  expansion  (dump)  combustor  work  carried  under  the  task  outlined  in 
Section  1  involves  interactions  with  personnel  in  the  Ramjet  Technology  Branch 
(RJT)  at  AFWAL/PORT,  Wright-Patterson  Air  Force  Base. 

In  addition.  Dr.  Edelman  chaired  the  JANNAF  workshop  entitled  "Turbulence  in 
Ramjet  Combustors"  held  in  conjunction  with  the  22nd  JANNAF  Combustion  Meeting, 
7-10  October,  1985  in  Pasadena,  California.  The  objective  of  this  workshop  was 
to  provide  a  link  between  current  basic  research  on  turbulence  to  problems  in 
ramjet  combustion  and  to  foster  the  utilization  of  relevant  information  on 
turbulence  in  the  design  of  ramjet  systems.  The  discussions  were  largely 
devoted  to  the  role  of  large  scale  structures  in  turbulent  mixing  and  combustion. 
Both  computer  simulations  and  experimental  work  were  presented  and  discussed  at 
the  workshop. 
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